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A gas of electrons in a one-dimensional periodic potential can be transported even in the ab¬ 
sence of a voltage bias if the potential is modulated slowly and periodically in time. Remarkably, 
the transferred charge per cycle is only sensitive to the topology of the path in parameter space. 
Although this so-called Thouless charge pump has first been proposed more than thirty years ago 
[1], it has not yet been realized. Here we report the first demonstration of topological Thouless 
pumping using ultracold atoms in a dynamically controlled optical super lattice. We observe a shift 
of the atomic cloud as a result of pumping and extract the topological invariance of the pumping 
process from this shift. We demonstrate the topological nature of the Thouless pump by varying the 
topology of the pumping path and verify that the topological pump indeed works in the quantum 
region by varying speed and temperature. 


Topology manifests itself in physics in a variety of 
ways [2-4] , with the integer quantum Hall effect (IQHE) 
being one of the best-known examples in condensed mat¬ 
ter systems. There, the Hall conductance of a two- 
dimensional electron gas is quantized very precisely in 
units of fundamental constants [5]. As discussed in the 
celebrated Thouless-Kohmoto-Nightingale-den Nijs pa¬ 
per [4], this quantized value is given by a topological 
invariant, the sum of the Chern numbers of the occu¬ 
pied energy bands. In 1983, Thouless considered a seem¬ 
ingly different phenomenon of quantum transport of an 
electron gas in an infinite one-dimensional periodic po¬ 
tential, driven in a periodic cycle [1]. This appears to 
be similar to the famous Archimedes screw [7], which 
pumps water via a rotating spiral tube. However, while 
the Archimedes screw follows classical physics and the 
pumped amount of water can be continuously changed 
by tilting the screw, the charge pumped by the Thou¬ 
less pump is a topological quantum number and not 
affected by a smooth change of parameters [1]. Inter¬ 
estingly, this quantization of pumped charge shares the 
same topological origin as the IQHE. The charge pumped 
per cycle can be expressed by the Chern number defined 
over a (1+1) dimensional periodic Brillouin zone formed 
by quasimomentum and time. Although several single 
electron pumping experiments have been implemented in 
nanoscale devices such as quantum dots with modulated 
gate voltages [8-10] or surface acoustic waves in order 
to create a potential periodic in time [11], the topologi¬ 
cal Thouless pump with sufficiently long spatial period¬ 
icity (which is crucial for achieving topological quantiza¬ 
tion [1]) has not been realized in electron systems. 

In this Letter, we report a realization of Thouless’ 
topological charge pump by exploiting the controllabil¬ 
ity of ultracold atoms in an optical superlattice. Dif¬ 
ferently from recent realizations of topological bands in 
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two (spatial or synthetic) dimensions [12-17] our exper¬ 
iment explores the topology of a (1+1) dimensional adi¬ 
abatic process, in which a dynamically controllable one¬ 
dimensional optical superlattice is implemented following 
the proposal of Ref. [18]. Topological pumping is seen as 
a shift of the center of mass (CM) of an atomic cloud mea¬ 
sured with in situ imaging. We extract the Chern number 
of the pumping procedure from the average shift of the 
CM per pumping cycle. The topological nature of the 
pump is revealed by the clear dependence on the topol¬ 
ogy of the pumping trajectories in parameter space. Our 
work introduces a new experimental platform to study 
topological quantum phenomena in adiabatic driven sys¬ 
tems. 

In our experiments, an ultracold Eermi gas of ytter¬ 
bium atoms ^^^Yb is prepared (see Methods A) and 
loaded into a dynamically controlled optical superlattice. 
Specifically, we construct a stationary 266nm-spacing op¬ 
tical lattice (short lattice) created by retro-reflection of 
a 532nm laser beam and a 532nm-spacing optical lattice 
(long lattice) with a time-dependent phase created by 
the interference of two 532nm laser beams, which is sta¬ 
bilized by a Michelson interferometer (see Methods B). 
As a result, these laser beams create the required [18] 
time-dependent one-dimensional optical superlattice of 
the form 

V{z,t) = -Vsit)cos^ (^) “ U(i)cos2 (ff - Ht)) > 

( 1 ) 

where d = 532nm is the lattice constant of the superlat¬ 
tice, Vs is the depth of the short lattice, Vp the depth 
of the long lattice, and 0 is the phase difference between 
the two lattices. In our experiments. Vs and Vp are con¬ 
trolled by the respective laser powers and 0 by changing 
the optical path difference between the two interfering 
beams with a piezo-transducer (PZT)-mounted mirror, 
which enables us to sweep 0 up to ^ IItt corresponding 
to more than ten pumping cycles. In the following, we 
use the lattice constant d as the unit of length and the re¬ 
coil energy Er = j(Smd?) as the unit of energy, where 
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FIG. 1. The Rice-Mele model, a. Schematic of the Rice- 
Mele model, b. A pumping cycle sketched (qualitatively) in 
(5-A space, c. Schematic of the continuous Rice-Mele (cRM) 
pumping sequence. The pink shaded packet indicates the 
wave function of a particular atom initially localized at the 
unit cell i. The wave function shifts to right as the pumping 
proceeds and moves the atom to unit cell z+1 after one pump¬ 
ing cycle. The blue dashed curve and the green arrow indicate 
the harmonic conhnement (not in scale) and an initial hole, 
respectively. 


m is the atomic mass of (see Methods B). 

We load ^^^Yb atoms into an array of one-dimensional 
optical super lattices, ensuring that they occupy the low¬ 
est energy band (see Supplementary Information S3), and 
slowly sweep (j) over time. The lattice potential returns 
to its initial configuration whenever (j) changes by tt, thus 
completing a pumping cycle. Since the lattice potential 
is periodic both in space and time, one can define energy 
bands and corresponding topological invariants such as 
the Chern number in a k-t Brillouin zone, where k is 
the quasimomentum and t is the time. We have ensured 
that the band gap never closes during the whole pumping 
procedure so ideally the atoms stay in the lowest band 
during the adiabatic pumping process. The phase sweep 
breaks time-reversal symmetry and the energy bands can 
acquire a non-zero Chern number. According to [1], the 
number of atoms transferred in one cycle is then propor¬ 
tional to the Chern number of the occupied band. 

The ability to tune all parameters of the lattice poten¬ 
tial (1) independently in a dynamic way offers the op¬ 
portunity to realize various pumping protocols. In the 
absence of the static short lattice, V{z,t) describes a 


FIG. 2. Observation of cRM pumping and sliding 
lattice pumping, a-b. In situ absorption images on the 
GGD before and after 10 cRM pumpings. c. One-dimensional 
optical densities (integrated along the x axis) before pumping 
(red circles, same data as a) and after 10 cRM pumping (blue 
diamonds, same data as b). d. The center of mass (GM) of 
the atomic cloud after up to ten pumping cycles. Red circles 
and blue open diamonds indicate the GM shift of the sliding 
lattice and the cRM pumping lattice, respectively. Error bars 
denote standard deviation of five independent measurements. 


simple sliding lattice which Thouless originally proposed 
[1]. Including the Vs term, one realizes a double-well 
lattice illustrated in Fig. 1. A pictorial understanding 
of this alternative pumping process is provided by the 
tight-binding Rice-Mele model [6, 20], 

^ - {J - S)albi^i + h.c. + A(dtdi - blbi 

i 

( 2 ) 

where and bi are fermionic annihilation operators in 
the two sublattices of the i-th unit cell, J ± d is the tun¬ 
nelling amplitude within and between unit cells, and A 
denotes a staggered on-site energy offset, as shown in 
Fig. la. We ignore the spin degree of freedom since we 
can neglect the interaction between the two spin compo¬ 
nents due to a very small 5-wave scattering length [21]. 

Figure Ic shows the schematics of our “continuous 
Rice-Mele” (cRM) pumping sequence. Sweeping the 
phase linearly in time as = TrtjT the hopping am¬ 
plitudes and on-site energies are modulated periodically. 

Our ab initio calculation shows that the cRM pumping 
scheme used in the experiment is topologically equivalent 
to the Rice-Mele model for atoms that reside in the lowest 
energy band since the Chern numbers are the same (see 
Supplementary Information S2). In the following, we will 
thus use the tight-binding Rice-Mele Hamiltonian to sim- 
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FIG. 3. Topological aspects of cRM pumping, a. Charge pumped during a simple cRM pumping (b), topologically- 
nontrivial pumping (c) topologically-trivial pumping (d), and negative sweep cRM pumping (e). The vertical error bars denote 
the standard deviation of three to five independent measurements and the horizontal error bar for the negative sweep data 
indicates a possible shift due to the nonlinearity and hysteresis of the PZT for negative sweep direction, b-e. Pumping sequences 
in the S-A plane (top), the corresponding band structures in the k-t Brillouin zone (middle), and the Berry curvatures of the 
pumping cycles (bottom). The indices w in the top figures indicate the winding number of each trajectory around the origin, 
f. Schematic pumping sequences of the cycles c and d. The phase (j) and the depth Vl of the long lattice are linearly swept 
alternatingly. The short lattice depth Vs is fixed to 30 Er and the pumping period T is 100 ms for both cycles. 


plify the discussion of the pumping sequence as a closed 
trajectory in the 5-A parameter plane (Fig. lb). 

Figure 2 shows the main results of our pumping ex¬ 
periments. Our stable absorption imaging system with a 
charge-coupled-device (CCD) camera enables us to accu¬ 
rately measure the shift of the CM of the atomic cloud 
after several pumping cycles (see Supplementary Infor¬ 
mation S4), as shown in Figs. 2a and b. The period T 
is fixed to 50 ms for the results shown in Fig. 2. One 
can clearly recognize the sizable CM shift along the 2 ;- 
direction. We plot the in situ CM positions of the atomic 
cloud after a few pumping cycles in Fig. 2d. The aver¬ 
aged CM shift per cycle {z{t) — z{0)) /{td) of the cRM 
pumping with (Vs^Vl) = {20,30)Eji is evaluated to be 
0.94(7) for t < 6T. This provides a direct measurement 
of the Chern number of the occupied energy band, which 
is consistent with the ideal value u = 1. As a com¬ 
parison, the observed average CM shift per cycle of a 
sliding lattice (Vs.Vl) = (0,40)F^i^ is 0.94(4), which is 
again close to the ideal value of = 1. Classically it is 
quite intuitive that the sliding lattice is able to trans¬ 
fer atoms because the potential minima are moving in 
space. However, even though the potential minima of 
the cRM pump {Vs^Vl) = (20,30)E'i^ are not moving 
in space as shown in Fig. Ic, the pumping is topolog¬ 


ically equivalent because of the same Chern number of 
the occupied band. The cRM lattice has the same abil¬ 
ity to transfer atoms residing in the lowest energy band, 
even though the pumping is achieved by a sequence of 
quantum tunneling events between the double-wells (see 
Supplementary Information S3). We attribute the satu¬ 
rating behavior of the cRM pumping for t > 6T to the 
effect of the harmonic confinement, whose variation can 
be comparable to the band gap for large CM shift [22] 
(see Supplementary Information S5). 

A striking feature of our pump is its topological na¬ 
ture. In particular, the pumped amount in the Rice-Mele 
model [7, 24] is directly related to the topology of the 
trajectory in the 6-A plane. It only depends solely on 
the winding number w of the trajectory that encloses the 
origin (5 = A = 0 (see Supplementary Information S2). 
Note that electron pumping in restricted nano-devices 
[8-11] is not topological since there the amount of the 
charge pumped per cycle instead depends on the area of 
the enclosed parameter space [25], which is the geom¬ 
etry but not the topology of the trajectory. To high¬ 
light the topological nature of Rice-Mele pumping, we 
investigate four distinct pumping sequences with trajec¬ 
tories shown schematically in Figs. 3b-e. In Fig. 3a, we 
plot the CM shifts of two cRM pumping schemes with 
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{Vs^yL) = (20,30)^i^ (Fig. 3b, e) and two amplitude- 
modified cRM pumping schemes (Fig. 3c, d). Evidently, 
the sequence which does not wind around the origin 
(Fig. 3d) results in no pumping, and those with wind¬ 
ing trajectories (Fig. 3b, c and e) result in finite pump¬ 
ing. This is direct evidence of the topological nature of 
the pump. Note that the band structure in the k-t space 
of the non-trivial pumping sequence (Fig. 3c) is identi¬ 
cal to that of the trivial pumping (Fig. 3d). However, 
the Berry curvature and the Chern number of the low¬ 
est band are different. This highlights the fact that the 
pumped charge is a topological quantity, which depends 
on the wave function but not on the band dispersions. 
Furthermore, we also performed the cRM pumping with 
a negative sweep of the phase (j){t) = —'KtjT^ which cor¬ 
responds to an opposite winding in (5-A plane and the 
cloud is pumped to the opposite direction even though 
the band dispersion remains identical to that of the for¬ 
ward sweep pumping (Fig. 3e). 


a b 



FIG. 4. Conditions for quantum pumping. The av¬ 
eraged CM shift per cycle after five cycles, a. Sweep speed 
dependence of the cRM pumping. The horizontal axis cor¬ 
responds to the period T of one pumping cycle. The inset 
shows instantaneous band maxima Ei of the first band and 
minima E 2 of the second band for {Vs,Vl) = (30,30)Fi?. 
The pumped amount saturates as the pumping speed slows 
down. The dotted curve shows an exponential fit with time 
constant (6.5 ± 1.2) ms. The error bars denote standard de¬ 
viation of three independent measurements, b. Finite tem¬ 
perature effect in the cRM pumping. The pumped amount 
approaches the ideal value as the temperature decreases The 
error bars denote the standard deviation of five independent 
measurements and the horizontal error bar indicates standard 
deviation of five temperature measurements before the lattice 
loading. 

A crucial requirement of the topological Thouless 
pump is adiabaticity, requiring that the band gap never 
closes during the pumping process and that the atoms al¬ 
ways remain in the lowest energy band. Figure 4a shows 
the pumping period dependence of the cRM pumping 
with the depths of (Vs, Vl) = (30, 30)F/^. The data sug¬ 
gest that a pumping time of T > 30ms is long enough to 
achieve a saturated pumped charge in this lattice poten¬ 
tial. This can be understood from considering Landau- 


Zener transitions to the higher band. The instanta¬ 
neous energy gap AE(t) = E 2 — Ei changes in time as 
shown in the inset of Fig. 4a. The probability of dia- 
batic transitions is given by the Landau-Zener formula 
P = p ^ (AE^ij2)‘^/h^{AE{t)). In our 

case, AFniin = l-OFi? and AE'max = ^OEr. Since the 
energy sweep speed is on the order of SAE'max/F, we 
find that 27rr ^ T/(6 ms), which is consistent with the 
observed result. 

We next check the temperature dependence of our 
pump. Figure 4b shows the pumped amount as a func¬ 
tion of the temperature of the gas before loading into 
the lattice of the depths (Vs,Vl) = (25,30)^^. The 
temperature is tuned by changing the sympathetic evap¬ 
orative cooling condition while keeping the same number 
of ^^^Yb atoms. One can see that pumping reaches the 
ideal value at the lowest temperature of 54(1 )nK. While 
the temperature given in Fig. 4b is not the temperature 
in the lattice but that before lattice loading, the observed 
temperature dependence exhibits similar behavior to that 
discussed in Ref. [18]: pumping efficiency reaches its ideal 
value once the temperature is lower than the gap energy. 

Having demonstrated topological Thouless pumping 
using flexible optical super lattice setup, the scheme can 
be extended to even more novel setups. For example, by 
choosing a special lattice laser frequency, one can create 
a spin dependent superlattice [26] which realizes Z 2 -spin 
pumping [27], a counterpart of Z 2 topological insulators. 
Another possibility is to change the ratio of the long lat¬ 
tice and short lattice wavelengths by tuning the angle 
of the interferometric lattice, which realizes a superlat¬ 
tice with incommensurate ratios and the Aubry-Andre 
model [28] with fractional pumping [29] and anomalous 
pumping [30]. Furthermore, introducing interaction ef¬ 
fect is feasible and will open a door for experimental 
exploration of the interplay of topological quantum phe¬ 
nomena and interaction and correlation effects. 

Note added. Recently, we became aware of similar work 
carried out by Lohse et al. [33] observing the topological 
Thouless pumping with bosonic atoms in Mott insulator 
state. 


I. METHODS 

A. A. Preparation of a degenerate Fermi gas of 

Since ^^^Yb atoms have a very short s-wave scatter¬ 
ing length of -0.15nm [21], we use sympathetic evapora¬ 
tive cooling with ^^^Yb atoms to obtain the degenerate 
Fermi gas of ^^^Yb [31]. After collecting Yb atoms in a 
magneto-optical trap using the intercombination transi¬ 
tion (556nm), the two isotopes are loaded into a crossed 
far-off-resonant trap (FORT) with 532nm light. Sympa¬ 
thetic evaporative cooling is performed by continuously 
decreasing the FORT trap depth. After blowing away 
the ^^^Yb atoms by using a 556nm laser which is res- 



















5 


onant only for we obtain a pure degenerate gas 

of ^^^Yb atoms with two hyperfine spin components of 
\F = \j2^mF = ±1/2) and a typical number of atoms of 
1 X 10^ for each spin. The initial temperature before lat¬ 
tice loading can be as low as 60 nK, which is ^ 50 % of 
the Fermi temperature in the FORT. The trap fre¬ 
quencies of the FORT at the end of the evaporation are 
I^TT = (145,20,161), where the x'- and y'- 
axes are tilted from the lattice axes {x and y) by 45°. 


B. B. Setup for the optical superlattice 


Our one-dimensional optical superlattice setup is a 
part of our optical Lieb lattice system [8]. To stabi¬ 
lize the phase (j) of the interfering 532nm-spacing opti¬ 
cal lattice (long lattice) we construct a Michelson inter¬ 
ferometer with frequency stabilized 507nm laser, whose 
optical paths are overlapped with the interfering lat¬ 
tice beams until separated from the lattice beams with 
dichroic mirrors before entering the experimental cham¬ 
ber. The 507nm laser beams are retro-reflected to form 
the interferometer just after the separation. We can con¬ 
trol the phase by tuning the PZT of the retro-reflection 
mirror of the z-axis 507nm laser, keeping stabilization of 
the optical path lengths of the lattice and 507nm beams 
via another PZT-mounted mirror in the common path 
as long as the phase sweep speed is smaller than 0.57r 
rad/ms. The short-term stability of the phase (j) is esti¬ 
mated to be 0.0077r. A phase drift of typically 0.057r per 
hour, is not a problem since the pumping only depends on 
the phase difference before and after pumping. The non¬ 
linearity of the PZT and the relative phase (p between the 
long lattice and the short lattice are calibrated via the 
matter wave interference pattern of a Bose-Einstein con¬ 
densate (BEC) of ^^^Yb atoms released from the super¬ 
lattice. The depths of optical lattices are also calibrated 
via pulsed lattice with the BEC of ^^^Yb. 


C. C. Calculations of the band structure and 
Chern number 

To predict the pumped charge of the experimental 
pumping protocols, we calculate the band structures and 
the Chern numbers of the one-dimensional Hamiltonian 
H{z,t) = —h V{z,t). The band structure is ob¬ 
tained by solving H{z,t) = F{k,t) in a 

plane wave basis. The Chern number is then calculated 
as 

2 pT n27rld 

z/= — / dt dkQ{k,t), (3) 

Jo Jo 

where ft{k,t) = dtA^ — dkAf is the Berry curvature and 
At{k) = is the Berry connection cal¬ 

culated using the periodic part of the Bloch wavefunction 

\ukit)) = IV’fe(i))- 
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Supplementary Information 

SI. POLARIZATION AND PUMPING IN A ID 
SYSTEM 

Thouless [1] calculated the current in a one- 
dimensional periodic potential using the Bloch wave func¬ 
tion and derived the formula for the quantized particle 
transport. Here we summarize the derivation of the same 
formula as the shift of a Wannier functions [2], which is 
more intuitive and relevant for the experimental realiza¬ 
tion of a finite and trapped system. 

Figure SI shows the schematics of our pump. We focus 
on a lattice site dX z = R. Using the full Bloch function 
of the lowest band IV^/c) = \uk)^ the Wannier function 

of the lowest band localized at the unit cell R is given by 


TKNN’s paper [4]. Therefore, the shift of the cloud after 
one pumping cycle is simply given as AP = where v 
is the first Chern number. In other words, the shift per 
cycle divided by d gives the Chern number. 

In the above derivation, we assume that the localized 
function can be written in the Wannier function consist¬ 
ing of the lowest band Bloch function. In the interacting 
case, a filled band is required for well-defined Wannier 
functions. However, in our experiments, we can ignore 
interactions and even if we have thermal excitations {i.e. 
holes), at the beginning of the pumping, each Wannier 
functions move independently and the center of mass of 
the cloud is robustly pumped as we observed. 


S2. TOPOLOGICAL EQUIVALENCE OF THE 
PUMPING PROTOCOLS 


1^) = ^ E = E 

^ k=—7T/d k=—n I d 

( 4 ) 

where N = L/d is the number of unit cells in the system, 
L the system length and d the lattice constant. The 
expected shift of the Wannier center from the lattice site 
P, or “polarization” [2] P at time t can be written as 


P(t) = {R{t)\z - R\R{t)) = 


= d 



dk , X 
—Ak{k,t). 


7r/d 

- y] {uk{t)\idk\uk{t)) 

— TV I d 


( 5 ) 


Here we used the relation dk\R) = 0. The quantity 
Ak{kR) = i {uk{t) \ dk \uk{t)) is called the Berry connec¬ 
tion and the integration on the right hand side of Eq. (5) 
is known as the Zak phase [3]. Thus, the spatial shift of 
the Wannier function from t = ti to t = t 2 or the change 
of the polarization is given by 


AP = 


' = P{t2) - P{ti) = d — [Ak{kR2) - Ak{k,h)] 

J-n/d 

(6) 

Using Stokes’s formula, this can be written as 


AP = -d 



dt{dkAt{kA) - dtAkikR)) (7) 


where At{kA) = i {uk{t) \ dt \uk{t)). Then the shift of the 
Wannier function after one pumping cycle is 


• nT nir/d 

AP = d— dt dk 

J 0 J—7r/d 






( 8 ) 

where = dk\uk). This formula is essentially the 

same as the one derived by Thouless [1]. Such an integral 
over two parameters in which the Hamiltonian is periodic 
should give an integer multiple of 27ri, as discussed in 


The simple sliding potential has Chern number u = 1 
for each energy band, therefore the pumped amount is 
proportional to the band occupation. This is a manifes¬ 
tation of the Galilean invariance of the lattice [5] . By lin¬ 
early tuning the lattice depth from (Vs, Vl) = (0,40)Ei? 
to the cRM pumping (Vs,Vl) = (20,30)E'i? with an in¬ 
terpolation parameter a G [0,1], the lowest two bands 
are adiabatically connected without gap closing as shown 
in Fig. S2a. The lowest two energy bands of the cRM 
pumping therefore also have Chern number u = 1. How¬ 
ever, a gap closes between the higher bands, and the 
third band of (Vs,Vl) = (20,30)E'i? lattice has Chern 
number z/ = — 1. The above discussion also shows that 
the topological equivalence of the cRM pumping adopted 
in the experiment and the two-band tight-binding Rice- 
Mele model [6] only holds for the lowest band, since the 
second energy band of the tight binding Rice-Mele model 
has Chern number z/ = — 1. 

For the Rice-Mele model, the pumping cycle defines 
a mapping from time to a closed path on the A — S 
plane. Focussing on the lowest energy band, we now 
show that the pumped amount is related to the wind¬ 
ing number of the path in the A — 6 plane [7]. In the 
periodic gauge the first term in Eq. (S7) vanishes, there¬ 
fore the change of the polarization is given by the path 
integral over the Zak phase, which is a smooth function 
except on the branch cut with a discontinuity of the size 
one. Eigure S2b clearly shows that the resulting change 
of the polarization will be finite only for paths that wind 
around A = (5 = 0. 

In the experimental setup, especially if Vl is compara¬ 
ble or larger than Vs, the modulated tunneling amplitude 
J A S changes drastically in the pumping sequence and 
J is not constant as in the Rice-Mele model. Moreover, 
hoppings beyond nearest-neighbors may be required to 
fully capture the band structures. However, these effects 
do not close the gap and thus do not change the topology 
of pumping sequence - a manifestation of the topological 
robustness of this system. Therefore, we can still charac¬ 
terize the pumping using 5 and A of the Rice-Mele model 
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FIG. S 1. Thouless pumping as a change of the polarization. 



FIG. S 2. Band gaps and Zak phases, a Band gaps of the interpolated pumping between {Vs,Vl) = (0,40)F^r and 
(20, 30)Er. AEij denotes the band gap between energy bands Ei and Ej. b Zak phase divided by 27r in the plane of A and S in 
the Rice-Mele model. The line of discontinuity can be chosen anywhere depending on the gauge phase choice of the eigenstate. 
The Zak phase or the polarization changes significantly when A changes the sign (See also discussion in S3). 


even in the experimental situations. 


S3. DETAILS OF 
DYNAMICALLY-CONTROLLED 
SUPERLATTICE SETUP 

In this section we explain the details of our 
dynamically-controlled optical superlattice setup and dis¬ 
cuss the band filling. The lattice laser configuration in 
the z-x plane is shown in Fig. S3a. The short lattice (a 
266nm-spacing optical lattice) in the ^-direction is cre¬ 
ated by retro-refiection of a 532nm laser (green arrow). 
The long lattice (a 532nm-spacing optical lattice with a 
time-dependent phase) is created by the interference of 
two 532 nm laser beams (yellow arrows), stabilized by a 
Michelson interferometer. As shown in Fig. S3a, the di¬ 
rection of the wavevector of this interferometric lattice or 
a “diagonal lattice” is not 2 ;-axis but {z — x) axis and the 
lattice has the lattice constant of d/ a/ 2 = 376 nm. More¬ 
over, the phase of the diagonal lattice moves to {z — x) 
direction as we sweep the PZT in the Michelson interfer¬ 
ometer. However, thanks to the x-direction confinement 


lattice Vio 64 created by a standing wave of a 1064 nm 
laser beam indicated by the red arrow in Fig. S3a, the 
hopping for the x-direction is suppressed and we can re¬ 
gard the 2 ;-direction projected component of this mov¬ 
ing diagonal lattice as a moving long lattice along the 
^-direction, which creates a one-dimensional superlattice 
potential for the 2 ;-direction (S3b). 

The setup contains another 532 nm laser to suppress 
the hopping for the ^-direction (H/) and realize an array 
of one-dimensional superlattice tubes (Fig. S3(c)). The 
depth of these vertical confinement lattices are Hfo 64 = 
3bEii and V]_ = IAQEr. Since we acquire a ^-axis inte¬ 
grated absorption image, we do not observe any differ¬ 
ence in pumping even if we lower U/ so that the tun¬ 
nelling time along the ^-direction is comparable to the 
pumping period. 

We also confirm the absence of the shift of the cloud 
along the direction perpendicular to the pumping direc¬ 
tion. For this purpose, we perform pumping experiments 
along the 2 ;-axis and the x-axis. In Fig. S4 we plot the 
CM position along the pumping direction ( 2 ;-axis for a 
and x-axis for b) and the perpendicular direction (x-axis 
for a and 2 ;-axis for b) for both superlattice configura- 
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FIG. S 3. Configuration of our superlattice setup.a. Laser configuration for the experimental realization of a 
one-dimensional optical superlattice. Black arrows indicate polarizations of the lattice beams, b. Lattice potential for 
(Fs, Vl, 1 ^ 1064 ) = (20,30,35) at 0 = 0 (upper) and 7r/2 (lower). The red dashed ovals indicate one-dimensional superlat¬ 
tice regions, c. Schematic view of our superlattice tubes array and the imaging axis. 

a b 




FIG. S 4. Absence of the shift of the cloud along the perpendicular direction, a The GM shift of the cRM 

pumping along the pumping direction (z, red circles) and the perpendicular direction (x, blue open diamonds) in the z-axis 
superlattice configuration (inset of a), b The GM shift of the cRM pumping along the pumping direction (x, blue open 
diamonds) and the perpendicular direction (z, red circles) in the x-axis superlattice configuration (inset of b). The lattice 
depths are (Fs, Vl) = (20, 30)Er for both a and b. 


tions (see insets in Figs. S4a and b). As one can see, 
there is no spatial shift along the perpendicular direc¬ 
tions, which is consistent with our system being an array 
of one-dimensional super lattices. 

To ensure that most of the atoms are loaded into the 
lowest band of the super lattice, we first ramp up the long 
lattice potential in 400 ms, whose band gap is much larger 
than that of the superlattice. We then ramp up the ad¬ 
ditional short lattice potential in 150 ms. The observed 
that at least 80 % of the atoms are mapped to the 1st 
Brillouin zone in our band mapping measurements. Con¬ 
sidering the non-adiabatic transition to the higher band 
during the band mapping, the actual occupation of the 
lowest band should be larger than this observed value. 

Moreover, we also confirmed that the atoms are always 
in the lowest band of the superlattice during one pumping 
sequence by using the “sublattice mapping technique” 


developed in our recent Lieb lattice experiment [8]. In 
this method, we first quickly change the lattice poten¬ 
tial to a “mapping lattice” with (Vs, V^qq^) = (40, 10)Er 
and a fixed relative phase 0 = 0 in 0.3 ms right after 
the pumping sequence. In this mapping lattice config¬ 
uration, the A-sublattice {a^} is energetically well sep¬ 
arated from the B-sublattice {bi} and the lowest two 
bands consist of the A- and B-sublattices, respectively. 
This maps sublattice occupations after the pumping se¬ 
quence to band occupations, which then can be mea¬ 
sured by band mapping techniques. Figure S5a, b and c 
show a schematic description of our sublattice mapping 
technique. In Fig. S5d, we show the sublattice mapping 
measurements for a cRM pumping sequence with the po¬ 
tential depths of (Vs,Vl) = (20,20)Er. As one can see 
in Fig. S5d, the band mapping result at t = 0, t = T and 
t = 2T are almost the same, namely the atoms in the 
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t=0 t=T/2 t=T t=3T/2 t=2T 

2fik,Kt 


FIG. S 5. Confirming occupancy of the lowest band. a. Sublattice mapping technique. If a wave function is localized 
at sublattice site the wave function is projected to the 2nd band of the mapping lattice and thus the occupation of the 
sublattice site ai can be measured via band mapping technique as the occupation of the 2nd band. The left image in a is a 
band mapped absorption image taken after 10 ms time-of-flight (five-times averaged), b. If the wave function is delocalized 
over both ai and bi, the wave function is mapped to both 1st and 2nd bands, c. If a wave function is localized at sublattice 
site bi, its occupation is mapped to the 1st Brillouin zone, d Sublattice mapping measurements (five times averaged). The 
horizontal scale of the TOF images are half of the TOF images in a-c. 


lowest band are still in the lowest band after pumping 
sequences. This result suggests that atoms localized on 
the B-sublattice {bi} can tunnel in the double-wells and 
then localize on the A-sublattice {ai}. The atoms are 
well described by the lowest-band Wannier functions of 
the superlattice. The result also suggests that the atoms 
are mainly localized at A-sublattice or B-sublattice in our 
deep-lattice cRM pumping sequence and that the tun¬ 
nelling only occurs when the Berry curvature has finite 
value. 


S4. DETAILS OF THE IN SITU IMAGING 

In this section, we discuss the stability of our in situ 
imaging. To evaluate a CM shift of an atomic cloud in 
real space, we measure the CM position from the in situ 
absorption imaging with a CCD camera (iXon, Andor). 
The CM position is derived from the Gaussian fitting 
to the integrated one-dimensional data. In our imag¬ 
ing setup, one pumping cycle, namely the cloud shift of 
d = 532 nm corresponds to 0.228 pixel on the CCD. The 
typical cloud width (FWHM) in the pumping direction 
is 3.7(5) pixels, corresponding to 16(2)d. 

Note that the CM shift after one pumping cycle is 
less than one pixel. However, since the Gaussian fit uses 
many data points (typically more than 20 pixels) as one 
can see in Fig.2c of the main text, we can reduce the 
fit uncertainty (la confidence) down to 0.3d for each im¬ 
age in Fig. 2 of the main text. Thus the error bars in 
the CM measurements mainly come from fluctuations of 



FIG. S 6. Stability of in situ imaging. The center of 
mass of the atomic cloud in the cRM pumping with depths 
of {Vs,Vl) = {20,26)Er. The error bars of the data in the 
region of 0 < t < T represent statistical errors of 30 to 32 
independent measurements and those in region of T < t < 2T 
are statistical errors of 9 to 12 independent measurements. 
The fitting uncertainty (~ 0.3d) is not taken into account in 
these error bars. 


the cloud position, probably due to the fluctuation of 
the initial atom loading. Although we cannot suppress 
this fluctuation, we can reduce the statistical error to be 
smaller than the fitting uncertainty by taking a lot of 
data. In Fig. S6, we show an example of such repeated 
measurements. The error bars in the region of 0 < t < T 
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FIG. S 7. Effect of the harmonic confinement, a Schematics of the Rice-Mele superlattice in harmonic confinement, b 
The local lattice potentials in three different regions (l)-(3) are shown in the upper part. The lower figures of b are schematic 
trajectories of the Rice-Mele pump in these three regions. 


represent statistical errors of 30 to 32 independent mea¬ 
surements and those in region of T < t < 2T statistical 
errors of 9 to 12 independent measurements. One can 
clearly see that the amount of shift is quite close to an 
ideal pump (zy = 1), indicated by a dashed line. 


S5. EFFECT OF HARMONIC CONFINEMENT 

In this section we discuss the effect of harmonic con¬ 
finement on the pumping. As illustrated in Fig. S7a, our 
system is finite and not uniform due to the harmonic con¬ 
finement. The spatial periodicity thus becomes worse as 
the CM of the cloud shifts away from the center of the 
trap. Moreover, since the local chemical potential also 
decreases as the CM of the cloud shifts from the center of 
the trap, the effective band gap becomes smaller and dia- 
batic transition occurs more easily - the edge of the cloud 
is more “metallic”. This can also be understood in the 


terms of the trajectory of the Rice-Mele pump as shown 
in Fig. S7b. The left (right) edge of the trajectory gets 
close to the degeneracy point as the cloud shifts to right 
(left), which suggest that the pumping should be satu¬ 
rated after a certain CM shift, depending on the strength 
of the confinement, while an ideal pumping is possible be¬ 
fore the saturation. In Fig. 2 and also Fig. 4 of the main 
text, one can indeed see that the pumped amount starts 
to saturate after several cycles. S8 compares the behavior 
of the saturation for two different confinement strengths, 
which clearly shows that larger confinement induces an 
earlier onset of saturation. 

This saturation point is also related to the value of 
the band gap. We actually observed that the saturation 
point becomes lower as we tune the long lattice depth Vl 
shallower corresponding to smaller minimum band gap. 
These observation suggests that the saturation point is 
determined by the competition between the band gap 
and the harmonic confinement. In other word, we can 
conclude that the atomic cloud are robustly pumped as 
long as the pumping system are protected by the gap. 
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FIG. S 8. Effect of the harmonic confinement in cRM 
pumping. The cRM pumping with depths of {Vs,Vl) = 
(20, 36)Er at different confinement strengths. In the case 
of stronger confinement pumping (entrap/S tt ~ 500 Hz, blue 
diamonds) starts to be saturated after a smaller number 
of cycles than that for the case of the weaker confinement 
(ct;trap/27r ~ 140 Hz, red circles). 




